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PREFACE 

This  research  was  sponsored  by  the  Advanced  Research  Projects  Agency 
under  the  monitorship  of  Mr.  William  Best  of  Air  Force  Office  of  Scientific 
Research.  Dr.  Shin-yi  Hsu  is  the  principal  investigator.  Research  scientists 
of  the  project  include  Dr.  Timothy  Masters  and  Ms.  Jane  Huang  of  Susquehanna 
Resources  and  Environment,  Inc. 

Mr.  Jack  Rachlin  and  his  associates  at  U.S.  Geological  urvey  served 
as  the  reviewer  of  the  effort.  It.  Colonel  James  Smith  of  AFOSR  served  as 
technical  adviser;  his  assistance  to  the  project  is  highly  appreciated. 

Phase  I  technical  report  has  been  read  by  Dr.  Mark  Settle  of  the  NASA 
Headquarters,  and  Dr.  Jack  Paris  of  Jet  Propulsion  Laboratory,  California 
Institute  of  Technology.  This  technical  report  has  taken  into  consideration 
their  comments  regarding  particularly  the  use  of  the  ordinary  geologic  maps 
as  the  ground  truth  information  of  the  LANDSAT  data,  and  background  infor- 
mation  on  textural  algorithms  for  lithologic  analysis. 
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•^Aimed  at  developing  image  processing  methods  for  rock  types  analysis 
with  LANDSAT  data,  numerous  experiments  were  conducted  using  supervised 
and  unsupervised  classification  techniques  under  the  general  concept  of 
texture  analysis  with  LANDSAT  digital  data  covering  two  geological  quads 

of  Nevada .  .  ,  .  ^ 

The  results  indicate  that  the  supervised  classification  method  is 
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, very  effective  in  the  extraction  of  granite  regions  when  (1)  data  were 
in  ratio  format,  (2)  feature  variables  included  both  tone  and  texture 
information,  and  (3)  the  classifier  is  capable  of  handling  non-normally 
distributed  data.  Classification  errors  occurred  when  there  exists  pixels 
of  non-granite  category  whose  spectral  and  textural  properties  are 
statistically  similar  to  that  of  granite  pixels.  Two  cases  of  errors  can 
be  noted:  Type  1  pixels  located  at  the  periphery  of  the  granite  regions, 
and  Type  2  pixels  located  far  away  from  the  core  of  the  granite  areas. 

To  reduce  the  error  rate,  an  unsupervised  classification  method 
based  on  the  concept  of  region  growing  and  texture  clustering  analysis 
was  employed  to  segment  the  scene  in  multiple  stages  and  thus  depict 
edge  patterns  by  the  scene  content  and  a  gradual  mathematical  generali¬ 
zation  process.  Identification  of  the  grantie  regions  becomes  a  labeling 
process  using  the  training  sets  information.  Since  the  Regions  algorithm 
is  based  on  an  additional  constraint  on  spatial  contiguity,  the  above- 
mentioned  two  types  of  errors  can  be  effectively  reduced  bacause  sharp 
edges  exist  between  the  granite  and  non-granite  pixels  in  the  study  area 
The  final  decision  regarding  the  delineation  of  the  granite  regions 
is  based  on  the  intersection  of  two  c 1  ass i f i ca t ion  maps  using  a  simple 
map  overlay  analysis.  The  result  yields  a  correct  c 1  ass i f i ca t i on  rate 
of  about  95  percent  based  on  a  visual  comparision  between  the  composite 
classif ication  map  and  the  ground  truth  information  given  in  the  U.S.G.S 
geological  map  of  the  study  area. 

To  improve  the  developed  techniques  for  lithological  analysis,  it 
is  recommended  that  additional  experiments  be  conducted  using  other 
regions  in  the  United  States  centering  around  the  following  tasks: 

(1)  developing  algorithms  for  merging  supervised  and  unsupervised 
classification  methods; 

I  (2)  finetuning  the  Region  algorithm  by  adding  subroutings  to  output 
digital  information  of  each  segmented  region; 

(3)  developing  a  color  prediction  model  for  rock  types  identification 
using  the  texture  and  tone  information  in  the  color  domain  with  a 
color  non i tor ;  and 

[b)  developing  change  detection  methods  for  monitoring  purposes  based 
on  the  extension  of  the  above  three  methods. 

The  above  discussions  apply  to  our  Phase  I  effort.  The  Phase  II 
i n ves tigation  is  designed  to  test  the  general izability  of  the  me t ho do  1 o 
gies  developed  from  the  Phase  I  experiments.  In  general,  it  has  been 
proven  that  they  are  indeed  genera  1 i zab 1 e  with  the  following  qualifica¬ 
tions: 

(l)  The  ratio  bands  are  not  an  absolute  requirement; 
j  (2)  Our  unsupervised  classification  method  has  been  improved  substantial¬ 
ly  to  the  point  that  it  can  be  used  as  a  smart  processor  for  extract¬ 
ing  alluvium  automatically;  and 

I  (3)  Quarternary  geologic  maps  are  more  appropriate  than  the  ordinary 
I  geologic  maps  for  serving  as  the  ground  truth  information  of  the 

LANDSAT  data. 

The  directions  of  future  research  should  be  centered  around  the 
'  development  of  smart  algor i thms  based  on  our  thorough  understanding  of 

the  physical  processes  by  which  the  terrain  units  were  derived.  ! 
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Executive  Summary 

It  has  been  determined  in  the  literature  on  seismology  and  geophysics 
that  the  recorded  seismic  wave  energy  from  nuclear  explosions  is  highly  depen¬ 
dent  upon  the  actual  yield  of  the  explosion  and  its  interaction  with  the  envir¬ 
onments  in  which  the  detonation  occurs.  These  environmental  factors  can  be 
characterized  by  the  depth  of  explosion  below  the  surface,  the  degree  of  coupl¬ 
ing  between  the  charge  and  the  adjacent  medium,  and  the  lithological  nature  of 
the  test  sites.  Therefore,  the  analysis  of  rock  type  at  the  test  sites  is  the 
first  step  in  nuclear  monitoring. 

The  LANDSAT  data  have  been  determined  effective  for  terrain  analysis. 

The  choice  of  the  LANDSAT  imagery  for  rock  types  analysis  at  the  nuclear  test 
sites  is  also  based  upon  the  fact  that  it  can  provide  world-wide  coverage  with 
repetitive  observatoins  for  monitoring  purposes.  And  in  certain  cases,  only  a 
combination  of  seismic  and  LANDSAT  imagery  can  yield  significant  information 
on  geology  and  tectronics  that  either  alone  would  not  provide  (Pavlin  and 
Langston,  1983).  The  goal  of  this  study  is  to  test  the  generalizability  of 
utilizing  LANDSAT's  digital,  multispectra L  information  for  rock  types  discrimi¬ 
nation  at  the  nuclear  test  sites,  based  on  the  texture-tone  analysis  algo¬ 
rithms  of  the  image  processing  systems  at  Susquehanna  Resources  and  Environ¬ 
ment,  Inc.  in  two  complementary  approaches:  supervised  classifaction  and 
unsupervised  training  methodology. 

The  experiments  were  based  on  two  subframes  of  LANDSAT  MSS  data  covering 
two  geological  quadrangles:  Site  1  located  in  the  Duffer  Peak  Quad,  Nevada, 
and  Site  2  uses  the  NE  quarter  of  the  Willow  Spring  and  Rosamond  Quads  of 
California.  In  addition  to  determining  the  general  capability  and  generaliza¬ 
bility  of  the  texture  analysis  algorithm  from  Site  1  to  Site  2,  this  study 
also  tests  the  appopriateness  of  using  a  general  geologic  map  as  the  ground 
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truth  for  judging  the  feature  extraction  capability  of  the  LANDSAT  data  based 
particularly  on  Site  2  where  a  quarternary  geologic  map  happens  to  be  avail¬ 
able  to  the  researchers. 

The  task  was  accomplished  by  using  two  separate  but  complementary  image 
processing  techniques.  The  first  technique,  a  supervised  classification,  was 
designed  to  extract  granite  regions  using  four  ratio  bands  (4/7,  4/6,  5/7,  and 
6/7)  for  the  Duffer  Peak  Site,  and  four  original  bands  for  the  Rosamond  Site 
based  upon  four  manually  selected,  but  automatically  pre-processed  training 
sets.  The  non-granite  regions  were  extracted  as  well  using  the  reject  cate¬ 
gory  of  the  classification  model.  The  second  method,  an  unsupervised  classifi¬ 
cation  procedure  based  on  the  concepts  of  the  stable  structure  of  scenes  (Hsu, 
1983),  was  designed  to  detect  the  contact  zones  between  granite  regions  and 
the  algorithm  using  one  ratio  band  (4/7  is  most  effective)  of  the  Duffer  Peak 
Site,  where  the  supervised  method  failed. 

For  Site  1  (Duffer  Peak,  Nevada),  the  final  granite  regions  were  defined, 
by  the  intersection  of  two  granite  images  produced  by  two  different  image  anal¬ 
ysis  techniques.  The  result  indicates  that  a  very  high  level  of  correct 
classification  rate — 95  percent  or  better — has  been  achieved,  based  on  an 
overlay  analysis  using  the  classification  result  against  the  geologic  map 
produced  by  the  U.S.  Geological  Survey. 

For  Site  2  (Rosamond,  California),  the  granite  regions  are  successfully 
extracted  using  the  raw  LANDSAT  MSS  information  based  upon  the  same  variables 
and  classifier  employed  in  Site  1;  and  thus  proving  the  generalizability  of 
the  developed  methodologies  for  rock  types  discrimination. 

On  the  appropriateness  of  using  the  general  geologic  map  as  the  ground 
truth  for  judging  the  performance  of  the  LANDSAT  system,  we  have  determined 
that  a  Quartenary  geologic  map  where  information  on  the  properties  of  surfi- 
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cial  (soil)  material  is  availabe  is  much  more  appropriate  than  the  general 
geologic  map  where  the  distributions  of  rocks  is  mapped  according  to  mainly 
the  interpretation  of  the  distribution  of  the  bedrocks  by  the  particular  field 
geologists.  Moreover  the  LANDSAT  system  is  really  not  designed  for  detecting 
sub-surface  material  although  the  sub-surface  information  can  sometimes  be 
inferred  from  the  surficial  expression  of  the  terrain  characteristics  as  cap¬ 
tured  in  the  LANDSAT  imagery  through  image  processing  and  analysis. 

Though  the  defined  task  of  extracting  granite  regions  has  been  success¬ 
fully  accomplished,  it  is  necessary  to  test  the  developed  image  processing  and 
analysis  techniques  using  additional  test  sites.  The  reasons  are  (1)  fine 
tuning  of  the  methods  are  usually  required  to  handle  diverse  patterns  of  litho 
logical  associations,  and  (2)  the  LANDSAT  imagery  can  be  exploited  further  for 
detecting  environmental  and  man-made  changes  before  and  after  nuclear  explo¬ 
sions,  but  it  has  not  been  fully  investigated  by  the  researschers  at  Susque¬ 
hanna  Resources  &  Environment,  Inc. 
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Applications  of  Texture  Analysis  for  Rock  Types  Discrimination:  Phase  II 

Section  A:  Introduction 

Ever  since  the  Soviet  Union's  detonation  of  its  first  nuclear  device  pro¬ 
totypes,  both  the  realities  of  an  arms  race  and  the  requirement  to  maintain 
scientific/technological  advantages  have  forced  the  United  States  to  expend 
significant  resources  in  monitoring  of  foreign  nuclear  tests.  Sophisticated 
technologies  that  have  evolved  about  the  framework  of  seismology  and  geophy¬ 
sics  have  made  significant  contributions  in  satisfying  the  national  require¬ 
ment  to  detect,  locate,  identify  and  yield-quantif y  world-wide  nuclear  detona¬ 
tions.  Yet,  there  is  room  for  improvement  using  non-seismic  methods,  particu¬ 
larly  in  the  area  of  yield  estimation.  To  this  end,  this  study  is  intended  to 

v  , 

develop  image  processing  and  analysis  methodologies  for  the  discrimination  and 
identification  of  rock  types  at  nuclear  test  sites.  The  rationale  of  this 
approach  is  based  on  the  fact  that  the  recorded  seismic  wave  energy  resulted 
from  nuclear  explosion  depends  on  the  following  environmental/lithological 
factors : 

(1)  the  actual  yield  of  the  explosion: 

(2)  depth  of  the  explosion  below  the  surface; 

(3)  the  degree  of  physcial  coupling  between  the  charge  and  the 
adjacent  medium;  and 

(4)  the  geological  nature  of  the  median  in  which  the  detonation  occurs. 
Indeed,  rock  types  analysis  is  the  first  step  in  yield  estimation. 

To  accomplish  the  goal  of  rock  types  discriminated  at  the  nuclear  test 
sites,  LANDSAT's  multispectral  data  were  used.  The  choice  of  the  LANDSAT 
imagery  is  based  on  the  fact  that  it  is  capable  of  providing  a  world-wide  and 
repetitive  coverages  and  thus  a  basis  for  monitoring  nuclear  test  activities. 
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The  thrust  of  this  study  is  to  exploit  the  digital  information  of  LANDSAT  data 
in  the  context  of  texture-tone  analysis  for  such  purposes. 

The  feasibility  of  the  SR&E's  image  processing  system  for  lithological 
analysis  has  been  proven  in  our  Phase  1  eftort  based  on  two  test  sites  in 
Nevada:  the  Antler  Peak  Quadrangle,  Nevada  at  the  scale  of  1:62,500  (ANA1 ) 
and  the  Duffer  Peak  Quadrangle,  Nevada  at  1:48,000  (ANA2)  as  analogs  to 
foreign  nuclear  test  sites.  Specifically,  the  first  site  (ANA1 )  was  used  as  a 
testbed  for  methodological  development;  whereas  the  second  site  (ANA2)  was 
designed  as  an  analog  area  for  extracting  granite  regions  used  in  the  Phase  1 
studies.  To  this  end,  another  site  was  selected  from  the  Willow  Springs  and 
Rosamond  Quadrangles,  California,  particularly  to  possess  these  properties  for 
testing  the  generalizability  of  the  developed  methologies: 

(1)  exposure  of  granite,  but  not  always  in  high  grounds; 

(2)  existence  of  contact  zones  between  granite  and  alluvium; 

(3)  availability  of  both  general  geologic  map  and  a  special  Quarternary 
geologic  map  where  information  on  surficial  material  is  given. 

Phase  II  effort  was  designed  to  test  the  generalizability  of  the  devel¬ 
oped  methodologies  used  in  Phase  I  studies.  To  classify  granite  versus  non¬ 
granite  regions,  two  complementary  image  analysis  techniques  were  employed. 

For  Site  1,  first  a  supervised  classification  analysis  was  conducted  to  delin¬ 
eate  granite  areas  based  on  manually  selected,  but  digitally  pre-processed 
training  sets,  and  to  reject  non-granite  regions  based  on  a  pre-set  statisti¬ 
cal  model/probability  level  for  identifying  pixels  which  are  significantly 
different  from  the  training  sets.  Second,  an  unsupervised  clustering  analysis 
based  on  SR&E's  Region  Growing  Texture  Clustering  algorithm  was  performed  to 
extract  granite  areas  by  region  growing  from  the  cores  of  the  granite  training 
sets.  The  final  definition  of  granite  regions  can  be  either  based  on  the 
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intersection  of  these  two  sets  of  "granite  maps,"  or  removing  alluvium  from 
the  supervised  classification  map  based  on  the  contact  zone  information  given 
by  the  unsupervised  classification  method.  For  Site  2,  the  analysis  is  much 
more  straightforward  since  a  very  accurate  classification  was  generated  by 
using  our  supervised  classification  methodologies  using  the  same  parameters/ 
variables  applied  to  Site  1. 


Section  B:  A  Brief  Review  of  Relevant  Literature  on  Lithologic  Analysis 

with  Image  Data 

Prior  to  1972  and  the  launch  of  LAND SAT,  pioneer  work  on  reflective  prop¬ 
erties  of  minerals  was  accomplished  by  Hunt  and  Salisbury  at  the  USAF  Cam¬ 
bridge  Research  Laboratories  (1970,  1973).  Their  study  and  explanation  of 
reflective/transmission  properties  of  both  minerals  and  rocks  in  the  visible 
and  near-infrared  regions  serves  as  a  basis  for  semi-automatic  rock  discrimina¬ 
tion  techniques  that  exploit  the  spectral  (tone)  parameters  of  multi-spectral 
imagery.  Rather  than  being  a  simple  empirical  result,  it  turns  out  that  miner¬ 
als  and  rocks  spectral  charcteristics  are  a  direct  function  of  the  physics  and 
theory  associated  with  crytal-field  theory  (Burns,  1970),  as  evidenced  from 
theoretical  and  laboratory  analyses  of  the  rocks  and  minerals  of  the  moon 
(McCord,  1968;  further  McCord,  et  al,  1972). 

Since  then,  scientists  at  the  U.S.  Geological  Survey,  Jet  Propulsion 
Laboratory  and  NASA/Goddard  Space  Flight  Center,  have  attempted  to  exploit 
LANDSAT  MSS  data  for  rock  *-ypes  analysis  as  evidenced  from  Goetz,  et  al 
(1973),  Goetz,  et  al  (1975),  Vincent,  et  al  (1975),  Rowan, et  al  (1976),  Rowan, 
et  al  0977),  Abrams,  et  al  (1977),  and  Podwysocki,  et  al  (1977).  Recent 
works  by  other  researchers  including  Lyon  (1977),  Lyon,  et  al  (1978),  Hunt 
(1977),  and  Siegrist  et  al  (1980),  also  emphasized  digital  processing  of  LAND- 
SAT  and  other  types  of  mult ispectral  scanner  data  for  optimal  combination  of 
spectral  channels  for  rock  discrimination. 

While  the  majority  of  the  work  cited  above  emphasized  rock  types  analysis 
and  identification  with  color  enhancement  techniques  with  LANDSAT  images,  our 
study  is  devoted  exlusively  to  extracting  rock  types  using  the  digital  informa¬ 
tion  of  the  LANDSAT  MSS  data  in  the  context  of  texture  analysis,  which  has 
been  largely  neglected  by  previous  researchers. 


Section  C.  A  Brief  Review  of  Texture  Measures  and  Testing  of  the  RADC/Hsu 

Texture-tone  Analysis  Algorithms 


tor  years,  texture  has  been  recognised  as  one  of  the  important  criteria 
for  identifying  objects  and  scenes  by  a  photointerpreter,  along  with  other 
variables  such  as  tone,  size,  shape,  associated  features,  etc.  Here  texture 
means  the  apparent  minute  pattern  of  detail  of  a  given  area,  described  ordi¬ 
narily  by  these  terms:  smooth,  fine,  rough,  coarse,  and  the  like.  In  digitaL 
data  processing,  texture  can  mean  the  spatial  distributions  of  tones  of  the 
pixels  of  a  given  area.  Its  attributes  have  to  be  specified  by  the  investiga¬ 
tor— a  specific  field  of  study  termed  texture  feature  extraction. 

Texture  analysis  is  a  rather  recent  but  rapidly  growing  field  of  inquiry, 
though  its  importance  relative  to  visual  perception  was  recognized  by  Gibson 
as  early  as  1950.  Over  the  past  20  years,  many  texture  measures  hve  been 
proposed.  This  body  of  literature  has  been  reviewed  by  Rosenfeld  in  1975.  in 
general,  these  measures  can  be  grouped  into  two  categories:  Fourier-based 
(power  spectrum)  features  and  statistical  features.  Furthermore,  it  has  been 

found  that  statistical  features  perform  much  better  than  the  others  (Rosen¬ 
feld,  1975). 


Haralick  (1975)  noted  further  that  there  have  been  six  basic  approaches: 
autocorrelation  functions  (Kaizer,  1955),  optical  transforms  (Lendaris  and 
Stanley,  1970),  digital  transforms  (Gramenopoulous ,  1973;  Hornung  and  Smith, 
1973;  Kirvida  and  Johnson,  1973),  edgeness  (Rosenfeld  and  Thurston,  1971)  and 
related  measures  (Schachter,  Lev,  Zucker,  and  Rosenfeld,  1977;  Lev,  Zucker, 
and  Rosenfeld,  1977),  structural  elements  (Matheron,  1967;  Serra,  1973),  and 
spatial  dependency  probabilities  (Haralick  et  al,  1973),  as  well  as  an  ex¬ 
tended  method  (Haralick,  1975).  In  general,  Weszka  and  Rosenfeld  (1975)  con- 
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eluded  that  statistical  features  perform  much  better  than  Fourier-based 
featu res • 

Texture  analysis  has  also  been  approached  from  the  human  perceptual  point 
of  view.  Indeed,  the  human  eyes  coupled  with  the  brain  are  very  effective  in 
identifying  and  interpreting  imagery  patterns.  The  only  drawback  is  the  slow¬ 
ness  of  data  processing  using  manual  operations.  Though  this  system  works 
empirically,  the  mechanism  by  which  visual  detection  and  recognition  is 
achieved  is  still  largely  unknown,  as  noted  by  Barlow,  Narsimhan,  and  Rosen- 
f eld  (1972)  and  Julesz  (1975).  This  field  of  study  has  been  termed  psycho- 
pietories,  as  summarized  in  Lipkin  and  Rosenfeld  (1970). 

With  respect  to  texture  perception,  Whitman  Richards  (of  MIT)  has  been 
conducting  experiments  under  the  sponsorship  of  the  Advanced  Research  Projects 
Agency  (1977).  He  has  concluded  that  most  uniform  textures  can  be  simulated 
by  three  or  four  variables,  provided  that  three  variables  contain  the  basic 
elemental  tokens  of  the  graphic  display.  His  approach  to  texture  perception  . 
has  employed  a  "generalized  colorimetry"  technique  analogous  to  that  used  so 
successfully  in  studying  human  color  vision.  Early  work  on  the  perception  of 
visual  texture  using  mainly  random  dots  includes  that  of  Pickett  (1967),  Polit 
(1976),  and  Purks  and  Richards  (1977).  Perceptually  based  texture  measures 
have  been  developed  by  Mitchell  et  al  (1977),  and  Hsu  (1977).  Thus,  we  see  a 
convergence  emerging  between  human  visual  processing  and  machine-oriented 
image  processing  methods.  As  will  be  discussed  below,  we  have  developed  a  tex¬ 
ture  analysis  system  that  indeed  is  capable  of  integrating  visual  processing 
into  a  machine-oriented  image  analysis  system. 

Under  the  sponsorship  of  U.S.  Air  Force/Rome  Air  Development  Center,  this 
author  developed  a  new  texture  measure  with  17  and  23  variables  derived  from 
(3  x  3,  Model  I)  and  (5  x  5,  Model  II)  windows,  respectively  (Hsu,  1977).  In 


the.  analysis,  the  window  moves  from  one  pixel  to  another  with  an  overlapping 
region  between  two  adjacent  pixels;  and  only  the  center  point  is  classified. 


TABLE  I 

The  Texture-Tone  Variables  of  Model  I 


Code 


Description  of  computational  funds 


1. 

MEAN 

2. 

STD 

3. 

SKEW 

4. 

KURT 

5. 

MDEVN 

6  • 

M  FT  CON 

7. 

MPTREL 

8. 

MINCON 

9. 

MINSQR 

10. 

M2NC0N 

11. 

M2NSQR 

12. 

MADAT1 

13. 

MAD AT 2 

14. 

MAD AT 3 

15. 

MBDAT1 

16. 

MB DAT 2 

17. 

MB DAT 3 

Average 

Standard  deviation 

Skewness 

Kurtosis 


) 


|  x_^  -  x|/n, 

|x  -  x  |/n, 
(x  -  x .  )/n 
lx.  -  x . I /n, 
-  x~! )  /  n 


where  x 


where  x 


i  and 


the  four  central  moments 

=  tone  value  of  individual  pixel 
=  mean 

=  tone  value  of  the  center  point 


j  are  adjacent  pixels 
i  and  k  are  second  neighbors 


(x. 

•xi  ■  vj"’ 

(x  -  X  )  In 

Mean  area  above  datum  1  (50) 
Mean  area  above  datum  2  (100) 
Mean  area  above  datum  3  (150) 
Mean  area  below  datum  1  (50) 
Mean  area  below  datum  2  (100) 
Mean  area  below  datum  3  (105) 


In  Model  I,  the  seventeen  texture  variables  are  as  follows:  (1)  through 
(4)  are  the  four  central  moments;  (5)  is  the  absolute  deviation  from  the  mean; 
(6)  is  the  contrast  of  the  center  point  from  its  neighbors;  (7)  is  the  mean 
brightness  of  the  center  point  relative  to  its  background;  (8)  is  the  contrast 
between  adjacent  neighbors;  (9)  is  the  sum  of  the  squared  values  of  (8);  (10) 

Is  the  contrast  between  the  second  neighbors;  (11)  is  the  sum  of  the  squared 

> 

values  of  (10);  and  (12)  through  (17)  are  the  mean  area  above  and  below  three 
datum  planes  having  tonal  values  of  50,  100,  and  150,  respectively,  on  a  scale 
of  0  for  black  and  255  for  white.  T'  a  code  names  and  computational  formulas 
of  these  17  variables  are  given  in  Table  1. 


t 


m  j 


I 


1 1 

•  i,  Model  11,  with  a  (3  x  ■})  design,  in  addition  to  the  above  17  vari- 
ibl.'s,  three  measures  are  extracted  to  characterize  the  oscillatory  nature  of 
t li,.  , ran  lines  obtained  along  both  the  x  and  y  axes  of  the  data  matrix;  thus, 
i -l  variables  are  available  for  analysis.  They  are:  (1)  sum  of  the  contrast 
values  from  peak  to  trough;  (2)  sum  of  the  distances  of  peak  positions  from 
the  origin;  and  (3)  sum  of  the  numbers  of  peaks  and  troughs  (see  Table  2). 
This  means  that  there  are  altogether  23  textun  variables  in  Model  II. 

TABLE  2 

Additional  Variables  in  Model  II 


Code 

Oescr i pt ion 

or  Formula 

18. 

XCONT 

(distances  from  peaks 

to  troughs)  along  x-axls 

19. 

XPKAK 

(peak  positions  from 

the  origin)  along  x-axis 

* 

20. 

XPANDT 

(number  of  peaks  and 

troughs)  along  x-axis 

21  . 

YCONT 

(distances  from  peaks 

to  troughs)  along  y-axis 

22. 

YPBAK 

(peak  positions  from 

the  origin)  along  y-axls 

• 

23. 

Yl’ANDT 

(number  of  pe  iks  and 

troughs)  along  y-axis 

In  our  work  for  the  Air  Force  Office  of  Scientific  Research  (Hsu,  1979 
and  1980;  Hsu  and  Burright,  1980),  we  developed  and  proved  the  existence  of  a 
perceptually-based  three-variable  texture  analysis  system  with  these  measure¬ 
ments:  (1)  average  tone,  (2)  the  first  neighbor  contrast,  and  (3)  mean  devia¬ 

tion  from  the  average  tone  in  our  Model  I  described  earlier. 


Section  t).  A  Brief  Note  on  the  Testing  of  the  RADC/llsu  Texture  AnaLysis 

Algorithms  for  Terrain  Analyses 

The  testing  of  the  effectiveness  of  the  RAI)C/Hsu  texture  analysis  algo¬ 


rithm  has  been  conducted  mainly  at  The  Pennsy Ivania  State  University  at 
College  Park,  Pa.,  and  at  the  State  University  of  New  York  at  Binghamton  in 
addition  to  the  current  research  performed  at  Susquehanna  Resources  and 
Environment,  Inc. 

The  Hsu  texture  measure  algorithm  was  implemented  at  the  Pennsylvania 
State  University  (Geophysics  Department)  by  G.  Pavlin  in  1979;  it  was  docu¬ 
mented  as  "Computer  Programs:  HSUDRIVK,  TEX  13  and  TEX  Several  theses 

have  been  'written  at  the  Pennsylvania  State  University  using  the  algorithm  to 
discriminate  lithologic  types  with  very  promising  results  in  arid  regions  (see 
Parker,  1980;  Ravenhnrst,  1980). 

At  the  State  University  of  New  York  at  Binghamton,  Kiracofe  (1983)  went  a 
step  further  to  tost  its  effectiveness  for  the  discrimination  of  vegetation 
types  in  addition  to  rock  types  in  the  Adirondack  region  of  New  York  State  and 
proved  the  algorithms  are  also  effective  in  terrain  analyses  in  both  arid  and 
temperate  regions. 

FinalLy,  we  should  note  that  the  U.S.  Air  Eorce/Rome  Air  Development 
Center  and  the  Defense  Mapping  Agency  have  implemented  these  algorithms  and 
used  them  for  both  terrain  analysis  and  target  cueing  for  years. 


this  next 
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Section  K :  Sunmviry  1 1 1  cl  Results  of  the  Phase  I  Kffort 

. 

I.  I  asks  Determined  for  the  Phase— 1  Kffort 

from  a  discussion  session  he Ld  among  Mr.  Best  of  AFOSR,  Col.  Lowrey  of 
DARPA,  Mr.  Raehl.in  and  his  colleagues  of  II. S.  Geological  Survey,  and  Dr.  Hsu, 
if  was  determined  that  the  tasks  of  Phase  I  effort  should  he  aimed  at 
answering  the  following  three  questions: 

1.  How  well  can  wo  map  the  granite  areas  versus  non-granite  regions 
using  our  supervised  and  unsuperv ised  classification  methods  in  the 
context  of  texture  analysis? 

2.  Hint  are  the  factors  affecting  the  classification  results — slope, 
drainage  pattern,  data  used,  methodologies  utilized? 

1.  What  an*  the  potential  contribution  of  image  processing  techniques 
and  methodologies  towards  the  discrimination  and  even  identification  of: 
rock  types  using  LAND SAT  data? 

For  data  analysis,  a  study  area  within  Duffer  Peak  Quadrangle,  Nevada  was 
selected  by  Mr.  Dempsie  of  U.S.  Geological  Survey.  Furthermore,  based  on  the 
geologic  map,  22  training  sets  were  selected  manually  to  cover  four  major  rock 
types:  (I)  granite,  (2)  tnetamorphic,  O)  volcanic,  and  (4)  unconsolidated. 

2 .  The  Data  Set 

To  remove  the  shadow  effect  of  the  original  LANDSAT  data,  and  to  extract 
information  from  four  MSS  hands  simultaneously,  the  following  data  sets  are 
gene  rated . 

(1)  First,  second  and  third  components  from  the  four  MSS  bands; 

(2)  Six  ratio  bands  from  the  four  MSS  bands:  4/5,  5/6,  6/7,  4/6,  4/7, 
and  5/7; 

C3)  The  first  component  map  from  4  seLected  ratio  bands. 

Therefore,  ten  derived  image  data  sets  are  available  for  analysis  in  addi- 
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tion  to  the  original  four  MSS  hands. 


The  location  of  the  training  sets  with 


respect  to  these  derived  data  sets  remain  the  same. 

5 •  image  Processing  and  Data  Analysis  Methodologies  Utilized 

To  analyze  the  relationship  between  the  selected  training  sets,  and  to 
c'  i«stfy  the  granite  areas  versus  non-granite  regions,  the  following  analyti¬ 
cal  techniques  are  utilized. 

;J  •  Extraction  ol  textnre-tone  information  of  the  training  sets  and  the 
entire  data  set . 

(sing  the  texture-tone  extraction  algorithm,  29  texture-tone-ratio 
variables  have  been  generated  for  any  given  pixel  from  four  multi-spectral 
bands  using  (3  x  3)  moving  grid.  They  are  composed  of  4  tone,  variables, 
l l  texture  variables  (1  from  each  hand),  6  ratio  variables,  and  6 
correlation  variables. 

For  data  analysis,  the  analyst  is  able  to  select  a  portion  of  the 
following  29  variables: 


BRIGHT* 

BRIGHT5 

BRIGIIT6 

BRIGHT7 

MINC0N4 

MINC0N5 

MINC0N6 

MINC0N7 

MDEVN4 

MDLVN5 

MDEVN6 

MDEVN7 

STDDEV4 

STDDEV5 

STDDEV6 

STD DEV 7 

MAXLIN 

LOG RAT 4 5 

LOGRAT46 

LOGRAT47 

LOG RAT 5 7 

LOG RAT 5 7 

LOGRAT67 

CORR45 

COR  114  6 

CORR47 

CORR56 

CORK 3 7 

CORR67 

Analysis  of  the  Training  Sets 

Based  upon  the  selected  variables  from  29-variable  system,  typically 
we  use  three  variables,  the  training  sets  will  be  analyzed  and  edited  so 
that  each  training  set  will  meet  the  following  two  criteria: 

(1)  single  mode;  if  two  modes  exist  in  one  training  set,  the  set 
will  be  split  into  subsets; 

(2)  extreme  outliers  arc  to  be  removed  based  on  a  statistical 
confidence  level. 


. .  . . . — .. _ 


c  .  J_) isc r i in i  nant  Analysis  of  the  Training  Sets 

Alter  the  training  sets  are  edited  or  preprocessed ,  they  will  be  ana- 

ly  nd  in  terms  of  how  close  they  are  between  pairs  of  training  sets  using 

the  means  of  selected  tone-texture  variables.  The  distance  is  generally 

2 

measured  by  statistical  distance  calling  Mahalanobis  D  with  or  without  a 

I og-determi nant  term. 

2 

While  the  D  distance  is  indicative  of  the  degree  of  similarity  and 
dissimilarity  between  two  training  sets,  the  analyst  usually  uses  a  con¬ 
fusion  matrix — classification  result  using  only  the  training  sets — to 
examine  how  well  these  training  sets  are  separated.  The  analyst  will  then 
decide  whether  he  should  proceed  with  a  classification  analysis  of  the 
entire  test  set.  in  general,  if  dissimilar  training  sets  are  confused,  a 
classification  analysis  should  not  be  conducted, 
d  .  .Supervised  Classification  Methods 

As  mentioned  earlier,  a  supervised  classification  analysis  can  be 
made  only  when  the  trailing  sets  are  well  separated.  To  achieve  this  goal, 
the  following  steps  can  be  taken: 

(1)  purify  the  training  sets  as  in  (b); 

(2)  change  the  location  of  the  training  sets; 

(3)  increase  the  power  of  the  feature  extractor  by  using 

(i)  more  texture-tone  variables,  and  (ii)  using  different 
spectral-band  combinations;  and 

(A)  increase  the  power  of  the  classifier  by  using  a  non-Gauseian 
model  if  the  data  are  essentially  non-mul ti-variate  normal. 

In  the  analysts,  we  have  done  all  these  image  processing  techniques 
’xcept  step  (2),  changing  the  Location  of  the  training  sets. 


The  Act  uni  Supervised  Classification  Procedures  Applied  Lo  the  Dul'fe 
Peak  Quad 

Step  I.  Four  ratio  bauds  (45 ,  5b ,  67 , 47  )  were  generated  as  the  basis 

Far  generating  29  tone-texture- lograti o-cur rein t ion 
variables  for  each  pixel. 

btep  2.  Preprocessing  of  the  manually-selected  Training  Sets. 

As  noted  training  sets  were  .selected  as  calibration  samples  of  four 
major  rock  types.  Using  information  from  bands  4/5  and  4/7,  these  origi¬ 
nal  22  sets  were  processed  into  25  sets  as  illustrated  in  page  13  of  the 
Phase  l  Report.  This  automated  preprocessing  technique  is  designed  to 
purify  the  Lraining  sets  by  Lwo  methods: 

L.  split  a  blmodal  distribution  into  two  sub-sets;  and 
7.  remove  outliers  of  the  distribution  from  the  training 
sets . 

Step  3.  Analysis  of  the  Preprocessed  Training  Sets. 

After  these  25  derived  training  sets  were  generated,  we  pro¬ 
ceeded  to  analyze  the  separation  pattern  among  them  using 
only  12  variables  from  the  29-variable  system.  These  12 
variables  are  composed  of  4  tones  and  8  texture  measures  from 
4  spectral  bands  as  follows: 

Tone  variables  are  from  the  4  spectral  bands; 

Texture  variables  are  the  1st  neighbor  contrast  and  the  Mean 
Deviation  measures  derived  from  a  (3x3)  moving  grid  from  all 
4  spectral  bands. 

The  result  of  a  confusion  analysis  was  given  in  Table  4  of 
the  Phase  f  Report . 


Step  4. 


Since  Ll  was 


evident;  f  fora  t  hin  anaLys  is  that  the  training 
sets  of  granite  were  confused  onLy  among  themselves  except 
Gl,  but  not  with  others  (see  Table  5  of  the  Phase  I  Report), 
we  then  proceeded  to  generate  only  training  sets  for  granite 
including  preprocessing  of  the  sets  using  G2,  C3,  C4  and  G5 
data,  (see  Table  3) 

This  step  is  to  set  up  a  supervised  classification  using  only  the 
training  sets  from  granite.  The  rest  of  rock  types  will  be  treated  as 
rejects  using  a  probability  cutoff  criterion. 

Step  5.  After  Step  4,  wc  proceeded  to  generate  a  classification  map 

using  only  7  feature  variables  in  the  classifier-4  tone  plus 
3  texture  measures  selected  from  the  original  28-variable 
system  as  follow: 

Brightness  4/5  1st  neighbor  contrast  5/6 

Brightness  5/6  Mean  Deviation  4/5 

Brightness  6/7  Standard  Deviation  6/7  (see  Table  4) 

Brightness  4/7 

SC(.p  6.  Printing  of  the  Decision  Map  as  given  in  Figure  I  of  the 

phase  I  and  Phase  1.1  Reports. 


TABLE  3 


Confusion  Matrix  Showing  only  the  Granite  Sets 


(Note  that  G2 

t  h rough 

G5  are 

confused  only  among  themselves) 

From 

Se  t  G 1 

02 

03 

04 

05 

Sum 

Cl 

34 

2 

3 

4 

3 

13.18 

0.00 

0.00 

0.00 

1.16 

G2 

0 

1.34 

0 

38 

71 

0.00 

48.20 

0.00 

13.67 

25.54 

87 

0  3 

0 

0 

61 

0 

1 .14 

0.00 

0.00 

32.80 

0.00 

61  .29 

94 

04 

0 

13 

1 

66 

56 

0.00 

.  30 

0.56 

37.08 

31.46 

76 

03 

0 

0 

L 

0 

1  .61 

0.00 

0 . 00 

0.62 

0.00 

99.38 

100 

TABLE  4 

The  7-Variable 

Texture 

-Tone  Analysis  System 

ST FP  1  VARIABLES:  BRI.GHT4  BRIGHTS  BRIGHT6  BRIGHT7  MINC0N5  MDEVN4  STODEV6 

Rank  of  e 

ach  set: 

TRANSTC2 

=  7  TKANSTG 3 

=  7  TRANSTG4 

=  7  TRANSTG5  =  7 

Confusion 

Matrix 

from  a 

Non-Gaussian 

Class 

tfier 

TRANSTG2 

TRANSTG3 

TRANSTG4 

TRANSTG5 

TRANS TO 2 

21 1 

0 

63 

5 

TRANSTG3 

0 

164 

2 

20 

TRANSTC4 

62 

2 

113 

2 

TRANSTG5 

12  ’ 

31 

19 

102 

Total  Correct  Classification  =  73.02  percent. 

Actual  Correct  Classification  =  73.02  percent. 

This  indicates  that  even  granite  sets  are  well  separated. 

f .  The  Experimental  Results 

From  three  experimental  data  sets,  it  can  he  concluded  that: 

(I)  For  rock  type  analysis,  data  Set  3  composed  of  these  4  ratio  bands: 


4/6,  A/7,  5/7,  6/7  is  most  effective.  Data  Set  1  with  the  4  original  4  LANDSAT 
MSS  bands  is  least  affective.  Figure  1,  a  decision  map,  indicates  that  the 
vast  majority  of  granite  areas  are  correctly  identified,  except 

(a)  tiie  granite  area,  within  which  the  training  set  G1  (which  was  not 
used  in  the  analysis  is  located,  is  largely  classified  ns  non¬ 
grin  ite,  and 

(b)  one  "metamorphi c  rock”  area  as  labeled  in  the  geologic  map  was 
largely  classified  as  "granite." 

These  two  regions  will,  be  investigated  further  using  our  unsupervised 
segmentation  algorithm  in  the  next  section. 

(2)  Regarding  the  classifiers,  our  non-Guassian  classifier  with  7  texture- 
tone  variables  is  superior  to  the  Guassian  classifier  no  matter  whether  it 
utilizes  7  or  16  text ure-t one  variables. 

(3)  There  is  little  difference  between  7-variable  Gaussian  classifier  and 
16-varlable  Gaussian  classifier  in  terms  of  the  confusion  matrix  using  the 
training  sets  data. 

(4)  In  terms  of  correct  classification  of  the  granite  versus  non-granite 
regions  (areal  distribution),  both  our  Gaussian  and  non-Gaussian  classifiers 
achieved  a  level  of  over  90  percent  hit-rate.  The  Non-Gaussian  Classifier  is 
slightly  better  than  the  Gaussian  Classifier  in  these  experiments. 

g .  Feature  Extraction  with  an  Unsnperv ised  Training  Approach 

( 1 )  Kxpcrimental  Design 

The  goal  of  these  analyses  is  to  extract  homogeneous  regions  in  the  study 
area  From  various  LANDSAT  ratio  bands  using  our  region-growing  texture  cluster¬ 
ing  analysis  algorithm.  Identification  of  the  granite  regions  becomes  a  label¬ 
ing  process  using  training  sets  information  and  other  related  statistical  and 


terrain  characteristics  data. 


It  was  our  intention  lo  iisp  the  results  from 
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this  unsupervised  classification  method  to  investigate  the  areas  of  misclassi- 
f  Lent  ion  by  tlie  supervised  classifier.  Seven  data  sets  were  used  in  the 

an  1 1 y s  is  . 

(?)  The  Results  of  the  Analyses 

By  examining  the  results  from  the  experiments ,  it  was  concluded  that  the 
results  from  Log  Ratio  of  Band  4  and  Band  7  (Figure  2  of  Phase  I  and  [I 
Reports)  yielded  the  best  overall  result.  Using  the  location  of  the  training 
sets  G 1  through  G5,  the  major  granite  areas  were  identified,  corresponding 
well  to  the  bedrock  regions  of  granite1.  Particularly,  by  comparing  the 
results  from  th Is  un, supervised  classification  against  the  supervised  classi¬ 
fication,  we  derived  that  (1)  the  rejected  granite  G1  area  can  be  delineated 
by  the  Region  algorithm,  and  (2)  the  confused  area  in  supervised  cLassifica 
Lion  between  G2  and  G3  can  be  discriminated  as  well.  Similar  to  the  super¬ 
vised  classifier,  the  Region  algorithm  faiLed  to  distinguish  the  bedrock  gran¬ 
ite  from  the  surficial  granite  in  the  are  near  training  set  G4.  As  it  was 
determined  in  other  experiments,  this  boundary  can  be  detected  in  the  analysis 
with  the  data  set  of  log  ratio  of  BamlA/Band  6. 

h.  Classification  Analysis  by  a  Combination  of  Supervised  and  Unsupervtsed 

Training  Approaches 

From  the  resuLts  given  earlier,  we  used  a  multiple  map  overlay  analysis 
to  delineate  the  final  granite  regions  ns  given  in  Figure  3  of  Phase  I  and 
phase  IT  Reports  with  the  following  conclusions: 

(1)  In  the  areas  where  training  sets  information  exists,  there  is  a  re¬ 
markable  correspondence  between  Figure  1  (supervised  classification)  and 
Figure  2  (unsupervised  classification); 

(2)  From  a  manual  editing  process,  we  can  place  the  granite  G1  area  Froin 


Figure  2  onto  Figure  1 ; 


(3)  The  areas  of  inlsclnssif  ication  in  Figure  1  — 

i.  region  between  G2  and  G3,  and 
ii.  pixels  located  outside  the  boundaries  of  labeled  granite 
regions  of  G1 ,  G2,  G3,  G4,  and  G5  in  the  northeast,  southeast  and 
southwest  quadrangles — 
can  be  removed  from  Figure  1. 

(4)  Since  there  is  no  information  regarding  ground  truth  in  the  area 
between  the  location  of  G1  and  G2,  we  used  the  result  as  given  in  Figure  1  for 
granite  identification;  and 

(5)  Comparing  the  results  as  described  above  in  reference  to  Figure  3,  it 
can  be  concluded  that  an  extremely  high  leveL  ot  correct  classification  of 
granite  and  non-granite  has  been  achieved,  it  should  he  noted  that  we  were 
able  to  edit  this  map  further  using  additional  ground  truth  information  and 
the  segmentation  results  given  by  other  ratio  bands. 
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Section  K:  The  Phase  II  Effort 

l.  Tasks  Defined. 

a.  To  test  genera l izahi iity  of  the  developed  methodologies. 

With  the  successful  results  obtained  through  the  Phase  I 
effort,  It  was  determined  that  the  tasks  for  the  Phase  II  effort 
should  be  centered  around  testing  the  generalizabllity  of  the 
developed  methodologies  for  rock  types  discrimination  using  another 
test  site.  To  this  end,  a  test  site  located  in  the  Willow  Spring 
and  Rosamond  Quadrangles,  California,  was  selected. 

During  the  Phase  I  effort,  both  Dr.  Smith  of  AFOSR  and  Dr.  Hsu 

felt  that  It  is  not  really  proper  to  use  geologic  maps  as  the  ground 

truth  to  evaluate  the  ability  of  the  LANDSAT  data  and  the  effective 

✓ 

ness  of  given  Linage  analysts  techniques  for  the  lithologic  analysis 
despite  the  Tact  that  it  was  employed  that  way  in  the  Phase  I  analy¬ 
sis  as  specified  by  the  researchers  at  U.S.  Geological  Survey  serv¬ 
ing  as  consultants  and  evaluator*  of  the  project.  Consequently,  a 
special  effort  was  made  to  select  a  site  for  which  both  the  standard 
U.S.G.S.  geologic  map  and  surficial  material  map  are  available  in 
addition  to  the  basic  criterion  for  conducting  a  generalizability 
study— comparability  between  the  first  test  site  and  the  second 
regarding  the  general  lithologic  characteristics  and  environmental 
conditions.  Based  upon  the  following  conditions,  the  Rosamond, 
California  Site  was  selected: 

(1)  Exposure  of  granite  at  a  much  lower  ground  for  testing  the 
effectiveness  of  the  developed  texture  analysis  algorithms 
on  extracting  granite  where  the . high-ground ,  well-exposed 
bedrock  condition  does  not  exist. 
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This  criterion  was  derived  front  a  general  feeLLng  that  tlu 
success  at  the  Duffer  Peak  Site  can  he  attributed  to  the 
condition  of  granite:  good  exposure  of  the  bedrock  due  to  a 
high-ground  location  lor  easier  detection  by  the  LANDSAT 
system. 

(2)  Availability  of  a  surficiaL  material,  map  for  the  test  site: 
Quarternary  geologic  maps  wore  located  for  several  sites  in 
the  II. S.  where  granite  can  be  found.  The  Rosamond  Site  was 
finally  determined  as  the  second  test  site  based  on  another 
criterion  on  the  similarity  in  climatic  conditions;  namely, 
they  are  Located  in  arid  climatic  zones, 
b.  To  evaluate  the  appropriateness  of  using  the  general  geologic 

map  as  ground  truth  for  LANDSAT  data  analysis. 

Although  wc  were  advised  by  the  U.S.G.S.  officials  serving  as 
advisors  and  evaluators  to  our  project  to  use  the  general  geologic 
maps  as  the  ground  truth  maps  for  the  study  areas,  both  Dr.  Smith 
and  Dr.  Hsu  could  not  totally  accept  this  approach  based  upon  the 
fact  that: 

(1)  The  general  geologic  maps  were  made  according  to  the  inter 
pretation  of  the  field  geoLogist(s)  regarding  the  distribu¬ 
tion  of  the  bedrocks;  but 

(2)  the  LANDSAT  system  was  designed  to  detect  primarily  the  sur- 
ficial  material,  and  thus  cannot  detect  sub-surface  material 
and  features  although  such  information  can  sometimes  be 
inferred  from  the  surficial  information. 

||(  e,  the  best  way  to  resolve  this  question  (dispute)  is  to 
treat  the  first  opinion  as  a  working  hypothesis,  and  test  it  against 


empiricaL  data.  With  the  availability  of  a  Quarternary  geologic  map 
for  the  second  test  site,  it  is  feasible  to  carry  out  this  task. 
c.  To  evaluate  the  contribution  of  the  unsupervised  classification 
method  to  tiie  extraction  of  lithologic  information. 

We  concluded  from  the  Phase  I  effort  that  both  supervised  and 
unsupervised  classification  methods  should  be  used  in  mapping  Litho¬ 
logic  features.  Since  then,  the  unsupervised  training  method  was 
developed  further  using  artificial  intelligence  approaches  by  the 
researchers  at  Susquehanna  Resource  &  Environment,  Inc.  According¬ 
ly,  a  portion  of  the  Duffer  Peak  test  site  was  used  to  demonstrate 
the  effectiveness  of  the  new  approach.  This  experiment  indeed  sup¬ 
ports  further  the  conclusion  derived  from  the  Phase  1  effort  that  a 
combination  of  supervised  and  unsupervised  classification  approaches 
are  needed  to  extract  lithologic  features  and  related  tectonic 
informal:  Ion . 

2.  The  Data  Sets 

A  frame  (256x256)  pixels  from  the  Los  Angeles,  California  landsat 
Frame  was  extracted  to  represent  the  second  test  site  located  at  Rosamond 
(Quadrangle),  California.  Since  Lt  is  determined  that  the  shadow  effect 
of  the  image  data  is  not  significant,  the  original  MSS  bands  were  used  as 
the  basis  for  data  analysis.  Note  that  ratio  bands  served  as  the  basis  of 
data  analyses  for  the  Duffer  Peak  Site  because  of  the  presence  of  shadow 
effect  In  the  image  data. 

-j.  Supervised  Classification  of  the  Second  Test  Site 

To  test  the  general inability  of  the  developed  methods,  identical 
procedures  used  in  the  phase  l  effort  will  be  employed  to  the  second  test- 


(a) 


,u.  The  fat  Lowing  .ectlon.  <le.ctll»d  the  methods  by  «Md>  «** 

types  were  extracted  from  the  Rosamond  Site. 

The  characteristics  of  the  training  sets 

As  mentioned  earlier,  we  treat  the  geologic  nap  oi  the  test 

site  as  the  hypothesized  ground  truth.  Hence,  the  training  sets 
were  selected  according  to  the  information  given  by  the  geologic 
map.  i gnre  is  a  portion  of  the  Willow  Spring  and  Rosamond 
Quadrangles,  California,  from  which  eleven  training  sets  were 
extracted  to  represent  five  rock  terrain  types  as  follows: 

(i)  A1  ,  A2  and  Ad  for  alluvium; 

(ii)  Cl,  G2,  C3  and  GA  tor  granite; 

( iii)  SI  for  schist, 

(iv)  T 1  and  T2  for  tuff;  and 

(v)  01  for  drybed. 

(b)  fioneratton  of  2’  cexture-tone  variables 

Fallow  tbs  some  procedure  used  in  the  Duffer  Peak  study;  twenty- 
nine  texture-tone  variables  were  generated  for  each  pixel  in  the 
(256x256)  frame  from  the  form  LANDSAT  MSS  bands  based  on  the  origi¬ 
nal  (3x3)  moving  grid  and  texture  analysis  model. 

(c)  Pre-processing  of  the  training  sets 

To  purify  the  training  sets  such  that  each  set  Is  uniform,  out¬ 
liers  in  each  training  set  (distribution)  were  removed  automatically 
using  Information  from  Bond  4  and  Band  7.  In  addition,  if  a  bimodal 
distribution  is  found.  It  will  be  split  into  two  subsets.  Analyses 
of  the  training  sets  Indicate  that  only  a  few  outliers  exist  In  each 

witnessed  from  the  following  data  (Table  >'• 


of  the  training  sets  as 


TABLE  5 

Pre-processing  of  the  Framing  Sets 


Training  Set  tD  No.  of  Pixels 


A 1 

256 

A  2 

2  56 

A3 

256 

Gl 

256 

G2 

256 

G3 

336 

G4 

256 

SI 

208 

T 1 

171 

T2 

209 

1)1 

(original)  No.  Pixels  after  Editing 

251 

246 

242 

249 

249 

324 

248 

242 

168 

209 


(d)  Selection  of  the  best  discriminators  from  the  29-variable 
system 

I)  determine  a  subset  from  a  large  number  of  variables  as  the 
host  discriminators,  a  stepwise  discriminant  analysis  is  generally 
employed.  However,  it  should  be  noted  that  the  results  of  the  anal 
ysis  are  highly  dependent  on  the  way  the  discriminant  function  is 
set  up.  For  instance,  in  one  case  the  best  discriminator  can  be 
designed  to  discriminate  the  closest  pair  of  the  training  sets;  on 
the  other,  it  can  be  devised  to  separate  training-set  classes  in¬ 
stead  of  individual  sets  in  the  classes. 

Since  the  goal  of  this  project  is  to  separate  granite  from 
alluvium,  the  step-wise  discriminant  analysis  utilized  here  follows 
the  second  approach— discrimination  of  classes.  The  analysis 
indicates  that  the  best  ten  variables  follows  closely  with  our 
original  concept  of  texture-tone  variables  as  follows: 

(1)  Bright  5  (Band  5  tone) 

(2)  Bright  4 

(3)  Bright  7 

(4)  Bright  6 


(5)  Standard  Deviation  fa 
(fa)  Mean  Deviation  5 

(7)  Lst  Neighbor  Contrast  7 

(8)  Mean  Deviation  4 

(9)  Standard  Deviation  ') 

(ID)  Mean  Deviation  7. 

(<*)  Confusion  analyses  of  the  training  sets 

To  determine  how  well  the  selected  variables  can  separate  the 
training  sets,  analyses  of  the  training  by  means  of  a  confusion 
matrix  is  generally  used.  This  is  achieved  by  classification  of 
individual  pixels  into  the  mean  vectors  oi  the  training  sets.  Table 
fa  is  the  confusion  matrix  from  the  10-variable  system  with  a  correct 
classification  rale  of  99.13  percent.  The  separation  between  train¬ 
ing  .sets  can  also  be  measured  by  means  of  the  Mahaianobis  distance 
(D“)  as  shown  in  TabLe  7. 

(f)  Classification  of  the  study  area  (test  site) 

To  be  consistent  to  the  procedures  used  in  the  Duffer  Peak 
experiment ,  we  used  the  identical  seven  feature  variables  to  analyze 
the  Lraiing  sets  again,  and  proceed  to  classify  the  entire  second 
test  site.  These  seven  variables  are: 

(L)  Bright  4,  (2)  Bright  3,  (3)  Bright  6,  (4)  Bright  7, 

(5)  lst  Neighbor  Contrast  3,  (6)  Mean  Deviation  4,  and 
(7)  Standard  Deviation  6. 

it  turns  out  that  these  seven  variables  is  a  subset  of  the  10- 
variable  system  derived  from  the  stepwise  discriminant  analysis. 
Corresponding  to  Tables  fa  and  7,  TabLe  8  and  Table  9  indicate  the 
confusion  matrix  ami  the  Mahalunobis  distance  between  training  sets, 


respectively.  Note  that  the  hLt-rate  with  respect  to  the  cla.ssiti 
cation  of  the  five  rock  types  is  99  percent. 

TABLE  6 

Confusion  Matrix  from  the  LO-Vartablc  System  (SR&E) 


A 1 
24  7 
0 
1  9 
0 
0 
0 
0 
0 
0 
0 
0 


A  2 
0 

246 

0 

0 

0 

I) 

0 

0 

0 

0 

0 


A3 

3 

0 

223 

0 

7 

1 

8 
0 
0 
0 
0 


G I 
0 
0 
0 

225 

18 

70 

30 

0 

0 

0 

0 


G2 

1 

0 

0 

5 

207 

0 

54 

0 

0 

0 

0 


G3 

0 

0 

1 

17 

0 

253 

1 

0 

0 

0 

0 


G4 

0 

0 

0 

1 

17 

0 

155 

0 

0 

0 

0 


D 1 
0 
0 
0 
0 
0 
0 
0 

242 

0 

0 

0 


SI 

0 

0 

0 

1 

0 

0 

0 

0 

204 

0 

0 


T 1 
0 
0 
3 
0 
0 
0 
0 
0 
1 

153 

8 


Classification  =  99.13  percent. 


TAB UK  7 


Mahalanobis  distances  from  Row  Set.  to  Column  Set  from  10-Variable  System  (SR&li) 


A 1 

A  2 

A3 

Cl 

02 

C3 

G4 

D 1 

SI 

Tl 

T2 

At 

0.0 

92.3 

10.6 

56.4 

29.3 

62.0 

24.9 

279.8 

195.7 

115.7 

170.0 

A2 

67.1 

0.0 

4  5.9 

53.  1 

50.3 

72.8 

44.7 

298.8 

168.0 

275.0 

470.6 

\3 

20.  1 

48.4 

0.0 

32.0 

17.6 

36.2 

13.1 

178.1 

122.  ! 

132.2 

252.0 

Cl 

48.  3 

44.2 

28.0 

0.0 

8.5 

3.6 

5.1 

537.1 

32.8 

34.6 

66.0 

S3.  7 

49.  i 

17.2 

1  i .  1 

0.0 

24.9 

2.  I 

469.5 

129.7 

149.5 

279.8 

0  1 

88.0 

57.9 

36.5 

3.5 

35.0 

0.0 

21.6 

629.5 

73.3 

44.7 

82.9 

04 

7^.5 

43.  1 

14.0 

14.8 

3.8 

30.4 

0.0 

415.7 

200.9 

173.1 

338.2 

ni 

4  31.2 

251.9 

435.6 

466. 3 

500.1 

4  37.7 

468.  7 

0.0 

412.3 

634.4 

775.  1 

SI 

202.7 

58.8 

110.6 

43.6 

88.8 

31.7 

74.9 

349.1 

0.0 

10.7 

22.2 

T! 

57.8 

136.7 

3  1  . 3 

33.4 

45.1 

36.3 

38.2 

908.9 

67.9 

0.0 

7.8 

T2 

53.2 

84.9 

36.9 

31.7 

32.1 

35.7 

29.9 

399.4 

37.7 

5.1 

0.0 

With  such  a  high  hit-rate  in  the  training  sets  classification, 
we  proceed  to  classify  the  entire  study  based  upon  these  (edited) 
training  sets.  The  result  is  given  in  Figure  4,  showing  the  deci¬ 


sion  on  the  distribution  of  granite  and  other  rock  types. 


TABLE  8 


Confusion  Matrix  of  the  Training  Sets  (Rosamond,  California  test  site) 

With  a  Seven-Variable  System 


Kent  ii  re  Variables:  BRIGIIT4  BRIGHTS  BR 1 GIIT6  BRTGHT7  MINC0N5  MDEVN4  STD  DRV  6 


Rank 

of 

each  set: 

A2 

=  7  A  3 

=  7 

G 1  = 

7 

G2 

=  7 

G3 

=  7 

1)1 

=  7  SI 

=  7 

Tl  = 

7 

T2 

=  7 

A 1 

A  2 

A3 

G  1 

G2 

G3 

G4 

m 

SI 

Tl 

A1 

243 

0 

6 

1 

0 

0 

1 

0 

0 

0 

A  2 

0 

246 

0 

0 

n 

0 

0 

0 

0 

0 

43 

19 

1 

0 

2 

0 

0 

0 

0 

0 

3 

G 1 

0 

0 

1 

225 

j 

17 

2 

0 

1 

0 

G2 

0 

0 

8 

26 

194 

0 

21 

0 

0 

0 

G3 

0 

0 

0 

82 

0 

239 

2 

0 

0 

l 

C.4 

0 

0 

7 

30 

50 

2 

159 

0 

0 

0 

D 1 

0 

0 

0 

0 

0 

0 

0 

242 

0 

0 

St 

0 

0 

0 

0 

0 

0 

0 

0 

203 

2 

Tt 

0 

0 

0 

0 

0 

0 

0 

0 

2 

149 

T2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11 

Classification  =  99  percent  with  respect  to  five  rock  types 


TABLE  9 


The  Mahal. tnobis  Distance  Between  the  Training  Sets  from  the  7-Variable  Systems 


Al 

A  2 

A3 

G 1 

G2 

G3 

C4 

01 

SI 

T 1 

T2 

Al 

0.0 

91.9 

10.1 

54.7 

27.0 

57.8 

22.8 

290.3 

195.6 

1  13.4 

161.1 

A  2 

60.6 

0.0 

42-4 

51.7 

44.9 

71.5 

40.4 

284.7 

168.4 

276.4 

467.1 

Al 

16.8 

46.9 

0.0 

31.9 

17.7 

35.6 

12.9 

169.7 

115.6 

113.2 

208.5 

01 

46.5 

4  3 . 0 

24.8 

0 . 0 

8.  1 

3.6 

5.3 

525.1 

33.2 

30.1 

60.2 

(!2 

61.9 

47.7 

15.9 

13.0 

0.0 

24.7 

2.1 

461.2 

128.4 

148.0 

272.4 

03 

86.6 

54  .8 

36.2 

5.0 

33.4 

0.0 

20 . 8 

602.8 

72.2 

38.  1 

70.2 

04 

71.0 

40. 1 

12.6 

14.6 

3.6 

29.3 

0.0 

395.7 

200.8 

174.1 

334.0 

D 1 

'*11.6 

240. 1 

412.8 

4  39.1 

478.0 

419.9 

448.9 

0.0 

383.1 

582.5 

717.6 

SI 

193.1 

56.6 

104.7 

38.4 

82.1 

29.2 

68.2 

341.7 

0.0 

9.8 

22.5 

T1 

49.5 

137.1 

28 . 2 

35.5 

38. 1 

37.4 

35.7 

832.6 

63.0 

0.0 

6, 5 

T2 

52.3 

84 . 7 

34 . 5 

31.6 

32.0 

35.1 

29.7 

394.9 

37.4 

5.2 

0.0 

4.  Discussions  on  the  Correspondence  Between  the  Decision  Map  from  the 

LANDSAT  Data  and  the  Genera L  Geologic  Map  and  the  Quarternary  Geologic 
Map,  Respectively. 

Since  the  goal  of  our  lithologic  analysis  is  to  separate  granite  from 
non-granite  regions,  we  will  concentrate  our  analyses  on  the  distribution 
of  granite  in  the  decision  map  against  that  delineated  in  (I)  the  geologic 
map,  and  (2)  the  Quarternary  geologic  map  as  follows, 
a.  A  comparison  against  the  1943  geologic  map 

To  provide  i  basis  for  comparative  analyses,  a  (10x10)  grid  was 


constructed  for  the  study  area.  Figure  5  is  a  portion  of  the  Willow 
Spring  and  Rosamond  Quadrangles,  California;  the  granite  region  was 
located  in  cells  along  rows  #4,  #5,  #6,  #7  and  # 8.  Alluvium  is 
located  mainly  on  the  northern  side  of  the  granite  region. 

Figure  6  is  the  decision  map  derived  from  the  LANDSAT  data 
based  on  the  above-discussed  training  sets  and  image  processing 
methods.  For  a  better  display;  it  is  broken  up  into  three  portions 
with  Figure  6a  showing  the  granite  region  at  the  middle  portion  of 
the  map,  and  Figure  6b  and  Figure  6c  indicating  the  upper  and  lower 
segments  of  the  classification  results,  respectively.  In  addition, 
the  outline  of  the  granite  regoin  from  Figure  5  was  traced  onto 
Figure  6  for  an  easier  visual  analysis. 

In  general,  the  computer  decision  map  (Figure  6a)  corresponds 
rather  wetl  to  the  geologic  map  In  terms  of  the  distribution  of 
granite.  Discrepancy  between  them  occurs  mainly  at  three  subregions 

(1)  Cells  #(6,2),  (6,3),  (6,4),  and  (7.4); 

(2)  Cells  #(4,7),  (4.8),  (4,9),  (4.10);  and 

(3)  Ceils  #(8,8),  (8,9),  and  (8,10). 

In  case  (1),  those  cells  were  designated  as  granite  in  the 
geologic  map;  whereas,  they  are  classified  essentially  as  alluvium. 

In  case  (2),  a  versed  pattern  of  case  (1)  is  witnessed. 

Case  (3)  is  similar  to  case  (1)  with  a  less  degree  of 
discrepancy . 

h.  A  comparison  against  the  Quarternary  geologic  map 

(l)  Note  on  characteristics  of  the  Quarternary  Geologic  Map. 

The  Quarternary  Geologic  Map  of  the  study  area  was  made  by 
researchers  (U.l.  Ponti ,  D.B.  Burke,  and  C.W.  Hedel)  of  U.S.G.S. 


33 

in  1981,  documented  ns  Open-File  Report  81—73  ,  and  entitled 
“Map  Showing  Quarternary  Geology  of  the  Central  Antelope  Vail  y 
and  Vicinity,  California."  Along  with  the  map,  "discussion"  ind 
"explanation"  sections  were  given  to  help  readers  understand  iow 
Lhe  map  was  compiled  from  field  work,  aerial  photos,  soil  sut 
veys,  and  other  source  material,  which  included  nineteen  reft  > 
uncus.  The  discussion  section  is  given  in  Appendix  1. 

The  researchers  noted  that:  "all  the  upper  Quarternary  map 
units  are  unconsolidated,  they  have  similar,  primarily  granitic., 
clast  lithologies,  and  they  remain  some  or  all  of  their  original 
depositions l  surfaces.  These  characteristics  distinguish  the 
deposits  irorn  older  Quarternary  and  pre-Quar ternary  formations 
of  diverse  lithology  which  are  weakly  to  firmly  consolidated  and 
deformed  and  which  preserve  none  of  their  original  deposttional 
surfaces . " 

From  this  explanation,  we  should  treat  upper-Quarternary 
deposi. t Iona  1  material  the  same  as  its  parent  material  in  the 
image  data  because  it  still  maintains  the  original  characteris 
tics  of  the  bedrock  lithology.  Since  eight  layers  of  Quarter¬ 
nary  datings  were  used  to  differentiate  the  relaLive  ages  of  the 
materiaL,  it  is  fairly  easy  for  the  users  to  understand  which 
map  units  belong  to  Upper  Quarternary,  and  which  is  Lower  Quar- 
Leruary .  In  addition,  each  map  unit  is  described  in  detail 
regarding  the  characteristics  of  stratification,  degree  of  con¬ 
solidation,  grain  size,  and  other  morphological  and  locational 
Information.  We  have  to  conclude  that  this  map  is  much  more 
appropriate  for  serving  as  the  ground  truth  information  Cor  the 

LANDSAT  imagery  of  the  study  area. 


(2)  A  comparison  between  the  LANDSAT  decision  map  and  the 
quar ternary  geologic  map. 

1'irst,  we  .si  ion 1  d  not1  that  the  Quarternary  geologic  map 
covers  only  a  portion  of  the  study  area;  fortunately,  the  major 
pint  ton  of  tin  granite  region  is  included  as  given  in  Figure 
7>  and  referenced  by  the  same  (10x10)  grid.  For  an  easier, 
comparative  analysis,  the  bedrock  granite  was  retraced  witli  a 
dark  line.  Now  let  us  compare  Figure  ha  against  Figure  7  in 
lemis  of  the  three  cases  of  discrepancies  noted  between  the 
computer  decision  and  the  general  (1943)  geologic  map  as  folLows 
!  n  Cells  #  (  ' ,  )  and  (6,3),  there  is  a  close  correspon¬ 
dence  between  these  two  maps.  This  cLearly  indicates 
that  the  (1943)  geologic  map  did  not  indicate  the  sur- 
ticial  material  of  that  locale,  which  is  composed  of 
sand  dunes  (Qds). 

Cells  #(6,3)  and  (7,3)  were  classified  as  alluvium. 

The  1931  Quarternary  geologic  map  shows  that  the  area 
has  been  highly  altered  by  human  actions  via  construc¬ 
tion  of  transportation  routes  such  as  highways  and 
railroads.  Therefore,  it  is  difficult  to  assess  that 
the  original  granitic  material  still  exists  today  and 
can  be  detected  by  LANDSAT. 

The  northern  border  of  granite  in  the  Quarternary 
geologic  map  along  row  4  is  located  further  north  as 
compared  to  the  1943  geologic  map  because  It  encom¬ 
passed  the  upper  Quarternary  deposits  (Q6m).  This 
boundary  actually  corresponds  well  with  the  decision 


map  of  the  LANDSAT  data. 
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Case  3.  The  material  Ln  the  upper  portion  of  Cel  Ls  //(8,8), 

(8,9),  and  (8,10)  were  classified  as  alluvium  by  the 
LANDSAT  data.  It  was  mapped  by  the  1943  geologic  map 
as  "granite."  In  the  Quaternary  geologic  map  it  was 
mapped  as  non-granite  bedrock  (or  outside  gr  m). 
However,  a  proper  Labeling  cannot  be  determined.  A 
field  cheek  is  needed  to  determine  the  precise  charac¬ 
teristics  of  the  material  at  this  location. 

5.  Conclusions 

From  the  analyses  given  Ln  previous  sections,  we  have  derived 
the  following  conclusions: 

a.  Our  developed  texture  analysis  algorithms  are  indeed  general 
izable  in  terms  of  their  effectiveness  in  extracting  granite 
versus  non-granite  lithologic  features  in  arid  regions, 
h.  Analyses  from  the  second  test  site  clearly  demonstrate  that 
Quarternary  geologic  maps  are  much  more  appropriate  than  the 
general  geologic  maps  for  serving  as  "ground  truth"  to  the 
LANDSAT  image  data.  Although  it  is  desirable  to  extract  bedrock 
information  from  the  LANDSAT  image  data,  it  is  definitely  not 
correct  to  conclude  that  the  LANDSAT  data  are  not  effective  when 
the  classification  map  (based  on  the  image  data  and  a  certain 
set  of  image  analysis  algorithms)  cannot  match  perfectly  to  the 
general  geologic  map  regarding  the  distribution  of  lithologic 
(bedrock)  features. 

c.  The  success  of  our  experiments  may  be.  attributed  to  three 
Image  analysis  nlorithms  that  are  not  generally  available  to 
other  researchers: 


/ 
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(1)  Training  set  pre-processor"  for  purifying  the  train¬ 
ing  set  data  so  that  each  set  is  uniform  internally. 

(2)  "Non-Gaussian  Classifier"  for  improving  the  classifies 
tion  result  when  the  training  sets  data  are  not  really 
normally  distributed. 

CD  "Region  Growing  Texture  Clustering  Algorithm"  for 
processing  the  image  data  in  an  unsupervised  training 
classification  mode.  This  algorithm  is  capable  of  detect¬ 
ing  the  contact  zone  between  granite  and  alluvium  regions, 
and  thus  extracting  the  feature  according  to  an  additional 
parameter  spatial  characteristics  of  the  interaction  among 
lithologic  units.  Although  this  algorithm  was  not  needed 
for  the  analyses  of  the  Rosamond,  California  test  site,  it 
was  necessary  for  editing  out  "granite  pixels"  that  are 
distributed  beyond  the  limit  of  the  granite  boundaries  in 
the  Duffer  Peak,  Nevada  test  site. 

With  newly  developed  capability  in  the  "Region"  algo¬ 
rithm,  we  will  demonstrate  its  ability  to  extract  alluvium 
Including  the  contact  zones. 
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Sect  Lon  G:  The  Duffer  Peak  Site  Revisited 

1 .  Introduction 

Tn  our  Phase  I  effort,  wc  employed  an  iterative  scene  segmentation  algo¬ 
rithm  to  depict  the  evolutionary  patterns  of  the  imagery  structures  of  the 
study  area,  starting  from  each  pixel  as  a  group,  and  ending  with  a  few  thou¬ 
sand  groups  in  the  (256x256=65,536)  frame.  Since  the  grouping  distances  are 
designed  to  progress  linearly  with  an  even  increment  of  one  unit  from  one  iter¬ 
ation  to  the  next:,  interior  pixels  are  grouped  first,  and  the  boundary  pixels 
(having  a  larger  neighbor  contrast)  should  remain  distinctive  for  a  while 
during  the  continuous  grouping  process,  revealing  the  contact  zone  between  two 
lithologic  units. 

To  a  certain  degree,  we  were  successful  in  identifying  granite  pixels 
(from  the  supervised  classification)  that  should  be  removed  from  the  original 
decision  map  because  they  were  located  beyond  the  "contact  zones."  Neverthe¬ 
less,  we  now  feel  that  a  better  algorithm  can  be  employed  to  identify  the 
contact  zone  directly,  instead  of  relying  on  a  continuous,  iterative  segmenta¬ 
tion  process,  the  original  region  algorithm.  The  new  algorithm  is  called 
"Edge"  algorithm  because  it  defines  feature  edge  (versus  spurious  edge)  accord¬ 
ing  to  certain  spatial  characteristics  corresponding  to  the  physical  proper¬ 
ties  of  modeled  features.  And  from  these  feature  edges,  we  should  he  able  to 
extract  the  alluvium  first  because  the  particles  of  alluvium  are  much  more 
uniform  than  those  of  granite,  and  the  contact  zone  between  granite  and  allu¬ 
vium  should  be  stronger  than  the  internal  edges  of  alluvium. 

2.  Demonstration  of  the  Capability  of  "Edge"  Algorithm. 

For  tills  experiment,  only  the  NW  quarter  of  Figure  1  (Decision  map  from 
the  supervised  classification)  will  be  used  in  view  of  the  fact  that  it  repre¬ 
sents  the  major  area  of  misclassif ication  by  the  supervised  method  due  to  the 
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. .  or  the  supervised  method  to  detect  the  cooler  tone.  Note  that  . 

supervised  classification  approach  i.  «««UU,  an  *•>*““  N°“ 

*  present  Figure  1  again  as  Figure  8  with  the  contact  tones  .ached  by  thichec 

llni„.  For  better  discussions,  the  NW  Quarter  is  subdivided  further  into  the 

£our  sections  as  indicated  in  Figure  8;  and  we  will  concentrate  our  analyses 

ln  the  lower  half  of  the  test  area,  particularly  the  SE  cell,  where  alluvium 

is  located. 

,et  un  recall  that  the  bedrock  boundary  and  the  granite  regions  nre 
narked  hy  "diagonal  strikes,"  rhe  no.erais  are  ciassfled  granite  pixel,  and 

"blanks"  are  non-granite  areas. 

„  ,  o,..  ceLL  of  the  NE  quarter  of  Figure  8,  It  is  essen- 

With  respect  to  the  Sh  ceu  oi  luk  i 

tinlly  a  non-granite  country  according  to  the  geologic  map;  however,  a  major 

portion  of  that  cell  was  classified  as  -granite"  according  to  the  LANDSAT  data 

„„d  our  algorithms.  Without  field  work,  it  is  difficult  to  seclude  which 

decision  is  right.  It  is  entirely  possible  that  they  ate  alluvia,  in  tecs  of 

the  grata  site  particles,  but  they  were  possibly  derived  fro.  the  parent  .ater- 

ial  (granite)  situated  at  higher  grounds. 

We  will  attempt  to  answer  the  above  question  particularly  by  means  of  a 
decision  »ap  generated  b,  o„r  "Edge"  algorithm  Particularly,  we  would  like 
to  point  out  that  the  decision  .ap  was  derived  b,  using  a  two-band  si.ulta- 
neous  segmentation  method!  the  first  band  is  represented  by  tbe  first  compo¬ 
nent  scores  of  tbe  original  four  HSS  bands,  and  the  second  band  is  the  ratio 

,  ,  a  RanH  7  The  result  is  shown  in  Figure  9  (a  region  map)  and 
between  Band  A  and  Band  /.  me  rcbun 

Figure  10  (an  edge  map  of  Figure  9). 

By  comparing  Figure  9  against  NW  Quad,  of  Figure  7  (classification  map 

with  a  supervised  classification  method),  we  can  immediately  notice  that  the 
contact  tone  between  the  granite  unit  and  the  alluvium  area  (as  indicated  in 
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L ho  geologic  map)  was  clearly  detected  by  the  Edge  Algorithm  with  a  very  high 
degree  of  accuracy.  In  addition,  the  alluvium  area  was  determined  as  one  uni¬ 
form  region,  whereas  the  granite  area  was  determined  as  a  composition  of  a 
great  number  of  highly  texturalized  local  features.  This  phenomenon  corre¬ 
sponds  closely  to  our  earlier  prediction  that  finer  particles  would  group 
together  sooner  than  coarser  particles. 

From  the  decision  map  of  Figure  9,  we  can  conclude  that  the  particles  in 
the  SE  quarter  of  Figure  9  and  its  continguous  part  in  the  NE  and  SW  quarter 
are  finer  than  those  in  the  rest  of  the  study  area  particularly  the  NW 
quarter.  However,  it  is  difficult  to  conclude  that  the  pixels  in  the  SE 
quarter  are  definitely  alluvium  in  the  general  sense  of  lithologic  classifi¬ 
cation.  From  the  terrain  data,  it  is  entirely  possible  that  those  fine  part¬ 
icles  were  derived  from  the  parent  material  called  granite.  If  those  mater¬ 
ials  are  of  upper  Quarternary  period,  they  could  still  maintain  the  property 
of  the  parent  rock.  Yet,  one  thing  is  sure:  the  contact  zone  between  granite 
and  non-granite  is  still  detectable  by  the  two-band  simultaneous  segmentation 
meth  d  of  our  "Edge"  algorithms.  This  contact  zone  is  clearly  shown  in  Figure 
10,  the  edge  version  of  Figure  9. 
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Section  H:  General  Conclusions 

From  Phase  I  and  Phase  II  research  efforts  on  rock  types  analysis  with 
LANDSAT  MSS  data  and  our  texture  analysis  algorithms,  we  have  come  a  long  way 
on  the  understanding  of  the  interaction  between  lithologic  material  with  the 
LANDSAT  sensing  system  on  one  level,  and  the  spatial  interaction  among  indi¬ 
vidual  pixels  and  groups  of  pixels  on  another.  Only  from  these  two  aspects  of 
"vertical-and-horizontal"  interactions,  would  we  have  a  better  handle  on  rock 
types  discrimination  and  identification  with  remote  sensing  technologies  and 
methodologies . 

At  the  beginning  of  this  project,  we  were  told  that  no  one  in  the  U.S.  in 
the  past  twenty  years  had  been  able  to  correctly  extract  granite  versus  non- 
granite  from  the  Duffer  Peak  test  site.  As  demonstrated  i.n  this  technical 
report,  researchers  at  SR&E  have  accomplished  the  task  that  has  been  extremely 
difficult  for  other  researchers  to  handle.  We  would  like  to  share  our 
thoughts  with  scientists  having  interests  in  lithologic  analyses  with  image 
data  on  what  we  have  learned  from  these  experiments  as  follows. 

1.  Technical  advisors  and  evaluators  to  our  project  have  insisted  on 
using  the  general  geologic  maps  as  "ground  truth"  information  for  the 
corresponding  image  data  because  in  many  cases  they  are  the  best  avail¬ 
able  material  to  the  research.  We  consider  this  position  as  too  rigid; 
and  properly,  geologic  maps  should  be  treated  as  one  of  the  many  informa¬ 
tion  sources  for  obtaining  real  ground  truth  that  can  be  detected  by  the 
sensors . 

2.  It  is  a  well-known  fact  that  the  characteristics  of  training  sets 
affect  significantly  the  final  classification  results.  However,  in  the 
past,  it  was  difficult  to  pre-process  the  training  set  automatically  with 
sophisticated  algorithms  because  of  hardware  limitations,  such  as  not 
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enough  core  capacity  to  perform  clustering  analyses.  Our  success  at  the 
Duffer  Peak  site  can  he  partially  attributed  to  good  training  sets  after 
they  were  pre-processed. 

3.  The  final  classification  results  can  also  be  affected  by  the  mathe¬ 
matical  model  that  the  classifier  uses.  If  the  training  sets  data  are 
Gaussian,  the  Gaussian  model  is  the  best;  if  not,  other  models  based  on 
non-parametric  distributions  are  better  than  a  Gaussian  model.  In  the 
past  and  even  today,  the  majority  of  researchers  are  still  using  Gaussian 
model  for  the  classifier  because  it  is  easy  to  program,  or  readily  avail¬ 
able.  however,  if  the  data  are  not  Gaussian  even  after  pre-processing, 
you  most  likely  will  increase  the  error  rate  by  five  to  ten  percent  in 
the  decision  map.  This  is  precisely  the  case  for  the  Duffer  Peak  site: 
Our  non-Gaussian  classifier  is  capable  ol  removing  a  substantial  amount 
of  "granite"  pixels  in  the  non-granite  area  that  are  present  in  the 
decision  map  generated  by  a  Gaussian  classifier. 

With  respect  to  the  Rosamond,  California  test  site,  both  the 
Gaussian  and  the  non-Gaussian  classifiers  produced  the  same  result. 

Thus,  either  model  is  acceptable. 

4.  Future  researcli  efforts  should  be  centered  around  developing  unsuper- 
vtsed  classification  methods  for  extracting  lithologic  features.  As 
demonstrated  in  this  project,  this  approach  can  be  effectively  used  in 
conjunction  with  a  well-designed  supervised  classification  to  determine 
the  contact  zones  among  various  lithologic  units,  and  thus  remove  the 
errors  made  by  the  classifiers  based  on  the  supervised  training  method. 

These  unsupervised  classification  methods  can  be  developed  into 
smart  algorithms  once  certain  decision  rules  are  implemented  into  the 
feature  extraction  processes.  Artificial  intelligence  can  therefore  play 
an  important  role  in  the  next  phase  or  generation  of  image  exploitation. 
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7 ni!,  map  it,  om:  nf  Lin  e  ■  in  the  U.5.  Geological  Survey  npen-f  i  les  thdl  emphas i ze  t  hi*  ndtui'O  anil  n ist rihut  inn  of 
upper  Quaternary  deposits  in  me  Antelope  Valley  ana  the  diljar.ent  canyons  ut  the  Iransverse  Ranges  ami  Tehar-hapi 
Mountains  111  south-central  California  (see  index  map),  the  area  covered  by  the  set  ol  three  maps  encompasses  about. 

J '>00  enuare  I  i  lome  tors  of  northern  Lus  Angeles  f.onlity  and  parts  of  San  Bernardino  60(1  hern  Counties,  excluding  Angeles 
a  0  r  r  °  dpits  of  towards  Air  force  Base,  lopoyraphy  in  the  area  ranges  from  rugged  sem-arid  mountains  and 
steep  canyon s^ to  broad  valleys  and  and  desert  flatlanos;  elovalinns  range  from  abon.  300  motors  in  val  oy  o,.;nds  ,o 
mire  tlidiH.ROU  meters  oo  mountain  peaks.  About  150,000  people  live  in  the  area,  with  most  of  Hit  popnlatlf.n  . 

-  are  mos, 

affected  by  land-use  decisions,  add  so  this  map  should  he  useful  to  planners  and  engineers  as  an  aid  in  assessing 
areas  subject  to  llash  floods,  f oimdat mu  and  drainage  problems,  severe  ground  motion  during  earthquakes,  and  other 
euolniiic  n.iyaids,  The  distribution,  age,  And  pattern  (if  faultinq  and  folding  of  the  depnsits  of  this  map  also  provide 
earth  sc  itii’  "  ‘  with  an  overvie*  of  t  lie  nature  of  sediment  deposition  and  deformation  in  nne  of  the  most  tectonically 

tic. live  reqi  of  the  world,  ,  .  ,  .  .  ,, 

This  I-  '•  desicjr.Od  js  si  regional  appraisal  of  th*  distribution  afi‘1  prooert  >•  s  of  late  Qu  si  ternary  materials.  It 
is  accurate  or  its  scale  and  purpose  as  an  aid  to  earthquake  liazaro  zonal  ion,  land-use  planning,  and  regional 
tectonic  analysis.  However*  it  snnulci  i> « *  considered  only  as  background  information  and  not  as  a  substitute  *oi 
largt-sca'e,  site-specific  studies  where  laud-use  and  engineering  decisions  ••■quire  more  detailed  geotechnical 

i nf onnat ion,  . ,  .  ..  r 

Uppir  (quaternary  diluvial,  colluvial,  lacustrine,  and  ^nlian  deposits  a  r«*  d  n  Teront  iated  on  the  mao.  Thus© 
materials  have  accuiau MU*<1  m  the  valleys  and  canyons  of  the  area  in  response  to  uplift  and  ernsjon  oj  the  Transverse 
RunytiS  and  TeMachapi  Mountains  ami  to  subsidence  of  the  Antelope  Valley  Basin  during  the  last .half-mill inn  years  01 
so.  Ail  tin*  upper  Quaternary  map  units  are  unconsolidated,  they  have  similar,  primarily  qranitic,  clast  lithologies, 
;mi  they  retain  some  or  all  nf  their  original  deposit iona I  surfaces.  These  character isi ic.s  distinguish  the  deposits 
from  older  Quaternary  and  pre-quaternary  formations  of  diverse  lithology  which  u«'"  weafcly  to  firmly  consolidated  and 
deformed  and  which  preserve  none  nf  their  original  depnsitional  surfaces. 

Aliuvidl  deposits  nf  Seven  major  episodes  of  deposition  are  the  most  widely  rxpnsed  upper  Quaternary  materials  m 
■the  did.  Cnnvletive  colluvium  with  generally  similar  textural  characteristics  and  alluvium  whose  texture  has  been 
modified  by  the  addition  of  windblown  sand  are  shown  on  the  map  with  distinguishing  patterns.  Materials  deposited 
during  the  high  stands  of  shallow  lakes,  alluvium  that  has  been  modified  by  the  addition  of  larqe  amounts  of  calcium 
•'urboiiate  around  the  Ijm  shoreline.,  and  dunes  of  uniform  sand  that  migrate  during  dry  lake  periods  occupy  the  valley 
.'ii. ns. 

...  determined  the  relative  ages  ot  the  upper  Quaternary  deposits  and  the  distribution  of  textural  facies  in 
deposits  by  compilation  ot  O.S.  Soil  Conservation  Service  soils  maps,  by  interpretation  of  aerial  photographs,  and  by 
Sludy  in  the  field.  Preliminary  maps  were  first  produced  by  compiling  soils  maps  of  Woodruff  and  others  (1  70), 

Using  their  descriptions  nl  the  major  soil  series  in  the  region,  we  were  able  to  identify  various  aip'S  ot  the  deposits 
and  tu  obtain  approximate  grain  size  distributions  in  soil  parent  materials.  We  could  dn  tins  because  for  deposits 
that  still  retain  some  ot  their  uriginal  surfaces,  the  degrees  of  profile  development  and  textures  of  the  soils 
directly  reflect  me  relative  ages  and  textures  of  the  deposits  upon  which  the  suils  formed.  Ihe  compilation  ot  sous 
mapping  tfien  served  as  i  guide  to  field  inspection  of  Soil  samples  from  channel  edge  exposures,  road  cuts, 
thousands  1  shallow  auger  holes  throughout  the  study  area.  From  this  information  we  determined  a  sequence  r*  sniis 
and  deposits  of  seven  distinct  ages.  Other  criteria  such  as  superposition  nf  deposits,  topographic  positin',  a"'! 
degree  nf  fan  surface  dissection  were  alsn  useful  for  relative  datinq  of  deposits,  particularly  m  areas  n  high 

relief,  structural  complexity,  and  windH>.  sand  veneers. 

1  no  grain  size  distributions  with  ->  Lin  cologic  units  differ  significantly  in  some  localities  from  the  oi  am  size 

interpreted  from  soils  descriptions.  aiK  K  re fore  relied  on  field  reconnaissance  and  the  examination  nf  seveia 

hundred  collected  and  sieved  samples  to  es-..  ri'  the  locations  of  t tie  facies.  Wide  variations  in  qrain  size  ovei 
small  distances  in  some  of  tin*  materials  make  arcifite  ‘K  1  ineat  ion  of  facies  Impossible  at  map  scale,  and  contacts 
between  facies  within  units  should  be  considered  as  1 ly  an"  • 

Radiometric  ages  of  the  uppei  Quaternary  units  arc  uiim 
the  area,  but  we  can  estimate  their  ages  from  stratigraphic  c  " 

Ihe  Pleistocene  Harold  Formation,  containing  land  mammal  fossus 
Q1  deposits  in  the  southeastern  part  of  the  area  (Noble,  1°F  ,  A. 
aye  of  the  oldest  Kanchoiahrean  fossils  is  460,000  years  (Repenning,  1980),  and  the  I 
deposits  on  this  map  is  thus  probably  no  older  than  about  400,000  yeirs.  The  ages  of 


te 1 >  local ed. 

mV  -•  1  le  material  is  very  i:i-  in  deposits  of 
„  j,  t. ..  i,.„  ,v  •  t  h  gated  deposits  el  sew  hero. 
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probably  no  older  tnao  about  400, UOH  yeirs.  me  ages  or  upper  Quaternary  mi".  *  " 
the  alluvia  I  sequence  are  estimated  on  the  assumption  that  the  depnsits  result  from  climatically  ^ntrollerl  fib's 
of  alluviation  that  are  essentially  synchronous  over  broad  regions  (Pont!,  1980;  Ponti  and  others,  1980).  1  be  uo it 
are  quite  extensive  ma  can  be  rocoqnized  alonq  both  the  Transverse  Ranges  and  Tehachapi  mounta  n  fronts  in  due  se 
tectonic  settings  and  across  the  various  microclimates  of  the  Antelope  Valley  and  adjacent  highlands-  Tmv  annea 

occurs 


itUIJIIM.  jail'hjj  O'lU  Uic  '.|.|...<U  . .  ...  .  ,  ..  -  ,  .  1  t  i.__  c 

have  their  origins  as  pulses  of  sediirvnt  produced  during  fluctuations  nf  climate  from  qlacial  in  interglacia 
and  cm  hr  tentatively  correlated  with  climatically  controlled  deposits  in  other  regions.  Good  correspondence  occur 
between  the  Antelope  Vall»y  units  and  the  Rivorbank  and  Modesto  formations  in  the  Ran  Joaquin  Valley.  01,  0?.  and  QJ 
deposits  appear  equivalent  to  tit*  upper,  middle  and  lower  members  of  the  RiveiTianl.  roimal  ion,  winch  liavi  OA  mat 
ages  from  45O.OO0  to  I30.0U&  years  (Marct.and  and  AUwardt.  1980).  Units  Q4,  Q6.  and  Q6  appear  equiva  eot  to  tin. 

Modes tu  Formation,  whicn  hal  an  age  estimated  to  span  the  period  frnm  90,000  tn  9,000  years  ago  ( Hare hand  and 
Allwaidt,  1980).  Lacustrine  deposits  (Opi)  and  calcium  carbonate  affected  alluvium  lQuca)  result  I rom  depos 
ana  groundwater  influences  of  pluvial  lakes  which  filled  the  Antelope  Valley  basin  during  the  most  recent  (P«t-03) 
glacial  periods,  units  q7,  present  stream  beds  (Qs),  and  sand  dunes  (Qds)  "re  in  part  historic  and  chanqe  each  seaso 
V,  1  tn  winter  raiiilall,  summer  t  lash  rlooos,  and  springtime  winds.  , 

contacts  between  units  beneath  upper  Quaternary  materials  are  compiled  and  simp  I  if  led  from  large-  ana 
intermediate-scale  mapping  by  tite  U.S.  Soil  Conservation  Service  (Woodruff  and  others,  1970),  and  from  Barrows  (1  7 

and  1980).  Harrows  ano  others  (197fi),  Boehy  (1977),  Dibblee  (I960,  1961,  1963,  1967),  dalins  and  Muehlherqer  (1954), 

K all  I e  (1577),  ami  Cable  and  others  (1976).  . 

boo  logic  structures  in  the  Antelope  Valley  are  from  unpublished  mapping  by  D.  (J.  Burke  and  t.W.  Iled"l,  those  io  the 
rift  zone  of  the  San  mamas  fault  are  from  recent  studies  by  the  California  Division  of  M  nes  and  Geo  oqy  B'l^ows, 

1 1/7/  aim  I960;  Barrows  and  others,  1976;  Derby  (  1977);  Kahle,  19/7;  Kahle  and  others,  1975).  those  in  tlu  ba.  lock 

l'g,3l-  "n(l  »h*  t,‘lb'' ’0Opy n avail abEf "to’ ^*1  iS '* ^oes  not r  (  ’’ 

permit  filly  1-gibb  iepicductiofl 


